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Abstract 

Results are presented for hadron spectroscopy with gauge groups SU{N) with N = 3, 5, 7. Cal- 
culations use the quenched approximation. Lattice spacings are matched using the static potential. 
Meson spectra show independence on N and vacuum-to-hadron matrix elements scale as y/N. The 
baryon spectrum shows the excitation levels of a rigid rotor. 



1 



I. INTRODUCTION 



Replacing the "3" of color SU(3) by "N"" and then taking N to infinity has a long history 
in the (continuum) phenomenology of the strong interactions, dating back to Refs. [l|, 0| • 

There is also a literature of lattice simulations applied to gauge theories with the group 
SU(N), for moderately large N. Most of it @-[3 



is directed at the properties of pure 
gauge theory I know of two papers on meson spectroscopy: Refs. (§, [Io| . They are done 
at smaller values of "large N," N = 2, 3, 4, 6. They are pretty standard lattice QCD 
spectroscopy calculations and reveal the N-independence of meson masses. Narayanan, 
Neuberger, and collaborators 11, 12| have measured masses and the pseudoscalar decay 



constant from simulations at much larger N (exploiting reduction, to simulate on smaller 
lattice sizes). 



But only Ref. [13j discusses large-N expectations for baryons, and it only makes compar- 
isons to actual lattice data for N = 3. So it seemed like an appropriate time to look at a 
lattice calculation of baryon spectroscopy at several values of N. 

Baryons in large N seem to be fascinating objects, either viewed as many-quark states 
[lij or as topological objects in effective theories of mesons 15, T^. There is an enormous 



(continuum) literature about spectroscopy and matrix elements for larg e A" baryons. As 



sorted early references include |l7H2l|. summarized in a review, Ref. 22] . Perhaps a lattice 
study might reveal something interesting? 



In fact, it does: general arguments [16|, [lj| state that the mass of an iV color baryon of 



angular momentum J should show a rotor spectrum: 

M(N, J) = NA + J ( J + ^ B (1) 

The formula applies to baryons made of an SU (2) isospin doublet of equal mass quarks. The 
parameters A and B depend on the quark mass and both should be some "typical hadron 
size," a few hundred MeV. The observation of Eq. [T] in simulation data is the main new 
result of this paper. Common expectations are that Eq. [1] is only true for small values of J, 
because then the terms have meaning as a good expansion in 1/N. However, I observe that 
it holds both for the bottom and the top of the multiplets where I did measurements. 

The first response of the lattice simulator to a proposal to look at baryons at large iV is 
undoubtedly negative: The cost of an iV-color simulation scales roughly like iV 3 , just from 
the cost of multiplying SU(N) matrices. Baryons are large objects so big volumes will be 
needed. Ordinary (N = 3) baryon correlators are noisy and because baryon masses scale 
as N, baryons in higher iV are probably noisier. And most of the literature involves some 
combination of small 1/N effects and exp ectation values of operators (£ ~ (B\0\B)), which 
are already hard enough for A^ = 3 [23] . Fortunately, large A" tests seem (mostly) not to 
require extrapolations to zero quark mass. 

Only for odd A" are baryons fermions, of course, so I am performing simulations of SU(N) 
gauge theories with AT = 3, 5, 7. 

I used the quenched approximation for all N. That means that, truly, I am not simulating 
QCD for any value of N. However, a first study does not quite demand the same level of 
quality as a later one might. I am interested in comparing simple observables in systems 
which differ only in N, but otherwise have the same UV cutoff (lattice spacing), lattice 
action, and physical volume. This is easy to set up in quenched approximation. 

Of course, the system which is simulated ought to have some connection to the real world. 
It is easy to justify a quenched approximation for large A" (although one could argue that 
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it remains an open question, whether such justification is correct). Quenched simulations 
for N = 3 are obsolete, but the effects of dynamical fermions on spectroscopy are really not 
all that large: away from the deep chiral limit they are small, order ten per cent for simple 
observables (compare Fig. 1 in Ref. {HI). 

The difficulty with using the quenched approximation comes when one wishes to convert 
a lattice number to MeV. The spectrum of quenched QCD is simply different from the 
spectrum of "real" QCD. I will compare dimensionless ratios of lattice quantities in my tests 
of scaling with N. Any conversions to MeV I make are only qualitative statements. 

One could take the position that one should test everything about large- N QCD at once, 
by doing simulations with dynamical fermions at several values of N. This just postpones 
the matching problem: the spectrum of the different N QCD's will be different. But perhaps 
that is not the right point of view. The real question one asks when comparing several values 
of N is whether dimensionless ratios of masses show some smooth behavior with N. And for 
that, one can begin simply, make observations, and then ask how they change as one does 
ever more realistic simulations. 

For all N's, I use lattices of size 16 3 x 32 points. I match the bare couplings so that the 
lattice spacings, as measured by various observables from the heavy quark potential, are the 
same, and then sweep in quark mass over similar ranges. This insures that finite volume and 
nonzero lattice spacing lattice artifacts will be similar across the board. I set the common 



scale using the shorter version of the Sommer |25| parameter, r\, defined in terms of the 
force F(r) so that r 2 F(r) = —1.0 at r = r x . The real-world value is r x = 0.31 fm 26], and 
with it my lattice spacings would be about 0.08 fm. 

Before going on, it's useful to set some definitions. The 't Hooft coupling is 

g 2 N = A (2) 



and the usual gauge coupling is thus 



The combination g 2 Cp, which appears in all perturbative calculations of renormalization 
factors, is equal to A(l — l/7V 2 )/2. For comparisons at fixed 't Hooft coupling, like the ones 
I present, this suggests that differences which might be perturbative (lattice-to-continuum 
matching factors, for example) will scale up to corrections of order 1/N 2 . 

The outline of the paper is as follows: The next section describes some technical problems 
I faced: using fat links, gauge fixing, and the construction of baryon operators. Only the 
third topic might be of interest to continuum physicists. Next, I present lattice results for the 
potential and for meson masses and simple matrix elements. The potential measurements 
show the extent to which lattice spacings are matched. The mesonic observables illustrate 
the N - independence of masses and the yN scaling of quantities like the pseudoscalar 
decay constant. Finally, in Section \W\ I show some results for baryon spectra. The major 
one is the presence of a rotor spectrum of excitations for two-flavor baryons. 
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II. TECHNICAL DETAILS 



A. Simulation techniques 



Mysimulations use a version of the publicly available package of the MILC collabora- 



tion 



27] . modified to generate gauge configurations and quark propagators at arbitrary N. 
Prior to this project, it had been extensively used in studies of N = 2, 3, 4 (28| . 

The gauge action is the usual Wilson action. Quenched simulations are performed using 
the standard mix of Brown - Woch microcanonical over-relaxation j29[ and Cabibbo - Mari- 
nari heat bath updates j30[, performed on all iV(iV — l)/2 SU(2) subgroups of the SU(N) 
link variables. It is known that simulations performed on SU(2) subgroups suffer critical 
slowing down at larger N, but the largest N is only 7 and this problem did not appear. 
Lattices are spaced 100 sweeps of the lattice apart, for later analysis. 

The spectroscopy was intended to be a typical SU(3) lattice QCD calculation writ large. 
This involved an improved fermion action and extended sources for hadronic correlation 
functions. To achieve these goals a number of small technical problems had to be solved. 
Many of them have been encountered in other large-iV studies, but perhaps my solutions 
are a bit different than what is found there, and might be worth reporting. 



B. SU(N) fat links; gauge fixing 

My lattice fermions are clover fermions with normalized hypercubic (nHYP-link) smeared 
links as their gauge connections 311]. The clover coefficient is fixed at its tree level value, 
c,sw — 1- This particular discretization is known to work well, with small scaling violations, 
in both ordinary QCD phenomenology and in beyond-Standard Model studies in SU(2), 
SU(3), SU(4). So the first technical problem is the construction of the nHYP link for 
arbitrary N. To describe the nHYP link in words, it is a local average of gauge connections 
over the hypercubes surrounding the link, which smears out the gauge field for the fermion. 
The specific problem to be solved is that the the fat link is the average of a set of paths 
which produces a sum of SU(N) matrices, call it E. The fat link is defined as 

v * = 7§i (4) 



This is the matrix which maximizes ReTrV^S^. The quantity Q x l 2 = l/v£t£ is com 



puted using the Cayley - Hamilton theorem as described in Appendix B of Ref. [32|. This 
construction involves finding the eigenvalues of Q, which is done using a Jacobi algorithm. 

As written, is an element of U(N), not SU(N). This is not an issue when it is used 
for the fermions since all that we care about is that our action should be gauge invariant, 
and under a gauge transformation both the thin and fat links transform the same way. 

State of the art spectroscopic calculations use extended sources as interpolating fields 
for hadrons. In this work configuration of link variables is gauge fixed to lattice Coulomb 
gauge and the source for the quark propagator is some spatially extended function. So, 
we need to gauge fix our lattices to lattice Coulomb gauge, rotating the links to maximize 
Si ReTr Ui(x), i labeling the spatial directions. This is done by finding a matrix V at 
each site of the lattice x which maximizes Re Tr VT) where E is the sum of forward going 
U^x)^ emanating from x minus the sum of TP(x — /i)'s terminating on site x. 
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In the author's SU(3) code, gauge fixing is done iteratively by sweeping through the 
lattice and, site by site, determining an optimal V by maximizing RelWEt. As in the 
case of configuration generation, the relaxation is done using the SU{2) subgroups of V. 
Unfortunately, for N > 4 relaxing on the subgroups suffers critical slowing down and it 
becomes impossible to carry out gauge fixing without performing an enormous number of 
iterations. 

This is a variation of a problem which has been previously observed (and solved) by 



several groups simulating large N gauge theories, with N > 10 |33l. I34j . There the problem 
is in the update step. Updating on 577(2) subgroups produces a simulation algorithm with 
a long autocorrelation time, which becomes longer with increasing N. The solution is to 



perform over-relaxation on the full SU(N) group, implementing an old idea of Creutz [35 



I do this a little differently than Refs. |33l . |34|. The maximizing V is given by Eq, HI which 



I have already dealt with while generating the fat links for the fermion action. The same 
application of the Cayley-Hamilton theorem gives V. As I already remarked, this V is an 
element of U(N), not SU(N). Now one needs SU(N) elements, so I must compute the 
determinant of V, extracting its phase <fi and performing an additional multiplication by the 
diagonal matrix exp(— 20/iV). Doing this makes gauge fixing no more expensive for SU(7) 
than for 577(3). 



C. Baryon operators 

Mesonic states are constructed in the usual way, by sandwiching fermion propagators 
with Dirac matrices. Baryons are a bit more complicated. In SU(3) it is common to use 
relativistic sources {e a bc[u a C^d b ]u c \ for example, for the proton). I don't know a nice way to 
generalize this to N > 3, so I built baryon states using non-relativistic quark model states. 
So far, I have only considered Nf = 2 flavors of quarks. In what follows, I will label the 
flavors as u and d, although one can give them different masses - in SU(3) the same operator 
can give the p or S, for example. 

In the MILC code, quark propagators are constructed in the Weyl (75 diagonal) basis, 
and are rotated to 70 basis. I keep the two "large" components of the propagator's source 
and sink spin indices to be the fields which contract the quark model states. There are two 
choices for projection (eigenvalues of 70 = ±1 allow forward going in time and backward 
going two-component quark propagators). One projection has the lightest positive parity 
state in a channel as the forward-going state, and the lightest negative-parity state as the 
backward-going one. The situation is reversed for the other 70 projector. In order to 
reduce noise, these two propagators are summed (actually subtracted), configuration by 
configuration. This produces a single correlator extending halfway across the lattice, which 
is a candidate for fitting to a single dying exponential. 

A generic iV— color baryon operator can be written as 

|5} = 6"^C {si} «...<|0). (5) 

The C's are an appropriate set of Clebsch- Gordon coefficients. The baryon propagator is 
a multiple sum over source and sink colors and spin coefficients of products of N quark 
propagators. This creates the possibility a large number of terms. Fortunately, many are 
redundant. 
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I condense states as follows. (I suspect that this is quite inefficient and that more efficiency 
is almost certainly possible.) I begin with states of definite u and d content. These states 
are eigenstates of J3, J, and J3. I then anticommute fermion fields into a "standard order," 
shown in Eq. moving from the left, spin- up u quarks, spin-down u quarks, spin- up d 
quarks, and spin-down d quarks. The epsilon symbol absorbs the resulting minus signs. 

Then, two examples are the I = J = Jz = N/2 state 

\B, I = J = N/2) = e abc - N u a t . . . uf |0> (6) 
and the SU(3) proton operator 

\p, t) = e abc ^(u a t u\dt - 14^4) |0> . (7) 
Wick's theorem says that the n— quark propagator itself is also a determinant: 

= (-)» J2 ^(P){ D ^^-y)D^{x-y)...D^(x- V )) 

P(l,2...n) 

(8) 

(P is a permutation of indices). The baryon propagator, then is a product of an up-quark 
determinant times a down quark determinant, summed over all the color combinations of 
the individual up and down quarks. Many terms are redundant in this product, as can be 
seen in an example, the SU (3) proton: The contributions from the e abc u^ui source terms are 
antisymmetric in both Eqs. [5] and El This means that it is only necessary to keep one color 
ordering (b > a, for example) for each pair of colors (a, b) in the baryon propagator, and 
each term can be reweighted by a factor of 2!. This generalizes straightforwardly so that a 
a term in any wave function with iV" spin-up u quarks, iV" spin-down u quarks, spin-up 
d quarks, Nf spin-down d quarks, picks up a restricted color sum, a single color ordering for 
each individual spin-flavor label and a multiplicity N¥\N±\N$\N?\. 

An extreme example of this pruning procedure is the the propagator for the / = N/2, J = 
N/2 state of Eq. El It is 

A(x,y) = (N\) 2 detM(x,y) (9) 

where M is the N x N matrix of propagators from a source color to a sink one. 

The lower the J, the more complicated are the wave functions. By SU(7), the J = 1/2 
states involve several hundred color combinations per spin configuration. This is getting 
quite unwieldy. (Since the nonrelativistic quark propagator is itself only a (2N) x (2N) 
matrix, there must be more redundancy!) But for now, I go on, naively. 

To summarize, flavor SU(2) wave functions and propagators for the various states are 

• I = J = N/2: as already described, there is one determinant of an iV x iV matrix, for 
a cost iV 3 . 

• Analogs of S* and H* states are 

\B, I = N/2 - 1, J = N/2) = y/(N-l)\(u$ . . . wf _1 sf |0) . (10) 

The baryon correlator is built of N 2 terms, one for each s color in the source and sink. 
Each is a determinant of an (N — 1) x (N — 1) matrix - the propagator of the up 
quarks. 
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I = J = N/2 — 1: We have an (N — 1) x (N — 1) matrix of u— quark propagators and 
iV + iV(iV — 1) color terms in the interpolating field for a cost of roughly iV 4+3 = N 7 : 



\B, I = J = N/2 - 1) = J • • • u ^ ld i «f~ 2 «f |0> (11) 

I = J = N/2 — 2: An (AT — 2) x (N — 2) matrix of propagators. The last term in the 
interpolating field has about A^ 3 color possibilities, for a cost of A^ 9 . 

\BJ=J = N/2 - 2) oc (m« . . . wf-^f-^f 



. . . TA-j. ^Ai, ^4. 
+u? . . . " 4 M f - 3 «f - 2 df ) |0> 



I = J = N/2 — 3: our lowest state if A" = 7; there are A^ 11 terms to evaluate. 
\B,I=J = N/2 - 3) oc («£... uf " 3 rff " 2 . . . df 

. "Z/^ G?^ G?^ 

~\~ 3ti^ . . . U-^ Xi^ It i G^ G^ 



-u a v .. <~ 6 <" 5 • • • <- 3 <" 2 • • • <) |0) 



(12) 



(13) 



Still to be constructed are three flavor (u, d, s) states. 

An obvious solution to the problem of increasing multiplicity is to prune states by doing 
an incomplete color sum. This will collide with another big problem: SU(N) baryon signals 
degrade as A" increases. For now, the only way to fight this is to collect larger data sets. 
Along the way, however, one can try to improve our signals by tactics such as averaging over 
the propagators with the same J and different m/s. (In practice, I combine the mj and 
—rrij propagators into a single correlator.) Another approach involves forcing a fit of several 
correlators which couple to the same states to a common mass. Ultimately, a variational 
calculation along the lines of what is done for excited state baryon spectroscopy might be 
necessary. 



III. RESULTS FOR PURE GAUGE AND MESONIC OBSERVABLES 

I present quenched results from three AT's, AT = 3, 5, 7. I have data sets of 80, 120, 
and 160 propagators, respectively The bare gauge coupling (3 is tuned to match potentials 
through the Sommer parameters tq or r±. (Recall r^F^o) = —1.65. t\ is the shorter-distance 
version of the Sommer parameter, t\F{t^) = —1.00. In the real world, ~ 0.31 fm.) 

In the present simulations, (3(N = 3) = 6.0175 (to match to previous work by Ref. [loj]). 
and then the N = 5 and 7 couplings were tuned to match r\. Couplings and derived 
quantities are recorded in Table IB The bare t Hooft couplings A turn out to be quite similar. 
Wilson-action simulations for A^ > 3 have a first order lattice-artifact transition, and I 
checked that my simulations are above it. Potentials for the three values of N are shown in 
Fig. [TJ They seem satisfactorily matched. 
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FIG. 1: Static potentials from our data sets: (a) STC7(3) at /3 = 6.0175, 16 3 x 32 volume; (b) 5*7(5) 
at /3 = 17.5, 16 3 x 32 volume; (c) SU(7) at /3 = 34.9. 12 3 x 32 volume. Effective mass fits for 
several values of t are overlaid. 



The combination r^^fa or r\\fa [a is the string tension) gives a dimensionless combination 
of the Sommer radius and the string tension a. The table shows that this quantity scales 
well. 

My spectroscopy is based on smeared Gaussian sources and p = point sinks. At each N 
I collected sets for several different values of the width of Rq for the source. These correlation 
functions are not variational since the source and sink are different. Thus, as Rq is varied, 
the effective mass (m e // is defined through fitting correlators at a single distance to a single 
exponential; with open boundary conditions for correlator C(i), m e // = log C(t)/C(t + 1)) 
can approach its asymptotic value from above or below. I observe that typically, as Rq rises, 
the mixed Gaussian - point correlators make their approach from above for smaller Rq, and 
from below for bigger Rq. An example of this behavior is shown in Fig. [2j (At large t, the 



S 



signal deteriorates - a characteristic feature quite familiar from many SU(3) studies. The 
noise is enhanced by the small data sets - 40 lattices - used to make the figure.) Then, 
rather than re-running the propagator code with yet another source, I can combine pairs of 
sources to produce a flat m e ff distribution using say 

C(t) = C(R = M) + fC(R = 8,t). (14) 

For SU(3), the optimal source is a linear combination of Ro = 4 and 6 sources, favoring 
larger R as the quark mass falls. For SU(5) I mix sources with R = 6 and 8. For SU (7) 
the Ro = 8 source produced flat effective mass plots across my mass range, and I did not do 
any source mixing. 

Tables HTllIVI contain the resulting spectroscopy. The values of the masses are highly 
correlated because they come from the same underlying configurations. Mass differences, 
which will be described below, are taken from jackknife averages of the data. 

Now I turn to results for mesonic observables. The critical hopping parameter k c is 
determined through the vanishing point for the Axial Ward Identity (AWI) quark mass, 
defined as 

dt MSfc t)° a ) = 2m, < P °( X ' *)° B > ■ ( 15 ) 

X X 

where the axial current A a ^ = ip^ '^75 (t° / '2)if) , the pseudoscalar density P a = ip r j^,{T a / '2)if) , 
and O a could be any source. Here it is my Gaussian shell model source. 

The mass is shifted from its free value in the usual way. To compare to usual expectations, 
recall the relation between g 2 Cp and A: we expect that the mass shift 

Sm=— 4 (16) 

2k c 

to show 1/iV 2 variation. This I roughly see; compare Fig. [3j Just for comparison over a 
wider N range, I show the older results of Ref. jsj, done with unimproved Wilson fermions. 

Next, we turn to spectroscopy, shown for mesons in Fig. [Hand baryons in Fig. |5j Data 
are plotted in terms of the AWI quark mass, in units of n, to make the x-axis the same for 
all N. The near independence of meson masses on N (when expressed in terms of a common 
variable) is apparent. Not so, for baryons! 

One is tempted to compare the chiral limit of the vector meson mass as a function of N. 
A simple linear fit to the data of Fig. Hlgives r^my = 1.58(2), 1.44(1) and 1.52(1) for N = 3, 
5, 7. With a nominal 1/n = 635 MeV, this is about 900 MeV to 1 GeV versus 770 MeV for 
the physical rho meson. One always has to be careful with quenched lattice results from a 
single volume and lattice spacing, but the number is completely sensible. It is close to the 
results of Refs. @, [lO 



(Using completely different methodology plus a different choice of physical input to set 
the lattice scale, the authors of Ref. [12] have a mass which is from thirty per cent to a 
factor of two higher than this. The difference between their result and the low-N ones will 
probably only be resolved by removing the many differences in methodology one at at time.) 

In what follows, I will replace the quark mass as the independent variable by the square 
of the ratio of the pseudoscalar to vector meson mass. The quark mass is scheme dependent. 
Of course, nothing I say is affected by this choice. 

Vacuum-to-meson matrix elements are expected to scale as y/N. I have looked at the 
pseudoscalar, vector meson, and axial vector meson decay constants. To be explicit, they 
are defined as follows: 

(0|u7o7 5 <i|7r) = m^U (17) 



9 



-I — I — I — I — I — I — I — I — r~ 



1.2 



1.0 -°> 



0.8 - 



0.6 



(a) 



$ ft 
m j I S I 



_l I I I I I I I I I I L_ 







5 10 
t 



15 



2.5 



2.0 



1.5 -a 



1.0 



-1 — I — I — I — I — I — I — I — I — I — I — I — r~ 



(b) 



- * $ - T 



_i i i i i i i i i_ 



= = [] 



5 10 15 

t 



. 1 1 1 1 1 1 


i 

(c) 






i 


i 


' 1 - 


~ $ 










[ 


I 


- T S 
■t*f Si] 


II BO 1 


) 
1 






( 


) 










- ( 


) 








t 


>< 


y < 


> 




■ , , , , 1 , 


1 - 




< 
f 




] 


, 1 ■ 



5 10 15 

t 



FIG. 2: Effective masses for the highest-spin baryon for different size sources, labeled by octagons 
for Rq = 4, diamonds for R = 6, and squares, R = 8. (a) 517(3), K = 0.1265; (b) 5*7(5) 12 3 
volume, k = 0.1275; (c) 5f7(7), 12 3 volume, k = 0.129. 

(so f n ~ 132 MeV) while the vector meson decay constant of state V is defined as 

(0\u7id\V) =m 2 v f v e l (18) 

and the axial vector meson decay constant of state A is 

(0\uj il5 d\A) = ttPjJau. (19) 

6j is a unit polarization vector. The lattice quantities f L are converted to continuum con- 
vention by f — f L (l — (3/c)/(4/c c )). I have left out the lattice-to-continuum Z-factor. With 
nHYP clover fermions, it is a few percent away from unity. 

The pseudoscalar decay constant is shown in Fig. |6j I show the dimensionless combination 
r ifps/y/N- It's n i ce to see the y/N scaling. Naive linear extrapolations give rif^/y/N = 
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FIG. 3: Additive mass shift versus 1/N 2 . My data are shown as octagons. The diamonds are the 
pure Wilson fermion data of Ref. jjj, just for comparison over a wider N range. 

0.154(2), 0.151(1) and 0.154(2) for N = 3,5, 7. The real world value is l/y/3 x 0.31 fm xl32 
MeV = 0.12, so the quenched decay constant at this lattice spacing is coming in about 15 
per cent high. Experts know that this kind of extrapolation is far too naive, to say nothing 
about comparing quenched QCD to the real world. Nevertheless, the answer is not absurd. 



Rcf. 24] shows a figure of with the ratio of quenched f % to its experimental result about ten 
per cent high. 

y/~N scaling for the pseudoscalar decay constant was first observed by the authors of 
Ref. [Tlj] . Their numerical result for the decay constant, translated to SU(3), is also high, 
but by 40 per cent. Detecting the origin of this difference will probably again require detailed 
numerical work. 

The same comparison is shown for f v , and f ai in Fig. [7J The a\ decay constant is quite 
noisy at small quark mass and I omit those results as untrustworthy. The y/~N scaling rule 
works well here, too. 



IV. RESULTS FOR BARYONS THE ROTOR SPECTRUM 

My baryon data in Fig. [5] shows masses which increase roughly linearly in N. All data 
shows that for the baryons, higher J does lie higher in energy. This is no surprise, so let's 
look deeper. Fig. [8] shows the splitting between the various members of each multiplet. It 
is extracted using a jackknife average of differences of the two baryon masses. 

I now demonstrate that the masses in Fig. [8] form a rotor spectrum. First, we can test 
the numerator of the rotor term of Eq. [1] N by N. This is done by looking at the ratio 
differences 

A(J M )= f'f-?,, (20) 
v i, 2, 3) M(N,Ji)-M(N,J 3 y v ; 
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FIG. 4: Meson spectroscopy in units of r\ from TV = 3, 5 and 7, plotted as squares, diamonds and 
octagons, respectively. Panel (a) shows the squared pseudoscalar mass. Panel(b) shows the vector 
meson mass. 
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FIG. 5: Baryon spectroscopy in units of n from N = 3, 5 and 7, plotted as squares, diamonds and 
octagons, respectively. For all N, higher J lies higher in mass: for N = 3 and 5 I have all states 
from J = 1/2 to J = N/2. For N = 7 only the J = 7/2, 5/2, and 3/2 states are shown. 
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FIG. 6: Pseudoscalar decay constant, matched from the different simulations using the r\ parameter 
from the potential and rescaled by 1/VN. Squares, N c = 3; diamonds, iV c = 5, octagons, N c = 7. 



0.20 




0.15 



0.10 



0.0 0.2 0.4 0.6 0.8 1.0 

\2 



( m ps/ m v)' 




0.05 

0.0 0.2 0.4 0.6 0.8 1.0 

v2 



( m ps/ m v)' 



FIG. 7: Vector meson and axial vector meson decay constants, rescaled by 1/v N. Squares, iV c = 3; 
diamonds, N c = 5, octagons, N c = 7. 
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FIG. 8: (a) Delta-proton (J = 3/2 — 1/2) mass splitting versus quark mass, in quenched SU(3). 

(b) SU(5) mass splittings: octagons for J = 5/2 — 1/2, squares for J = 5/2 — 3/2, and crosses, for 
J = 3/2 - 1/2. In large-N, Eq. [TJsays that the SU(5) J = 5/2 - 3/2 splitting is supposed to be 
equal to the SU(3) J = 3/2 — 1/2 mass splitting, and also to the SU(7) J = 7/2 — 5/2 splitting. 

(c) SU(7) mass splittings: octagons for J = 7/2 — 3/2, squares for J = 7/2 — 5/2, and crosses, for 
J = 5/2 - 3/2. 



for which the constants (A, B) cancel. The result is shown in Fig. O I just plot one 
mass difference as a function of the other one and compare the data to a straight line of 
zero intercept whose slope is given by Eq. [201 The rotor spectrum is confirmed for all the 
members of the N = 5 and 7 multiplets I observe. 

Second, one can look at the J = 3/2 — 1/2 splitting, and check the N dependence for the 
bottom of the multiplets: 

M(N, 3/2) = M(N, 1/2) = ^ (21) 
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FIG. 9: Mass differences in the SU(5) and SU(7) multiplets, panels (a) and (b) respectively. For 
SU(5) I show a line whose slope is the ratio A(5/2, 3/2, 1/2) (recall Eq. [20]) . The line has slope 
3/8. For SU(7) I show a line whose slope is A(7/2, 5/2, 3/2) = 5/12. 

Given the states I have recorded, this can only be done for N = 3 and 5. Rescaling the mass 
difference by iV/3 exposes B. The result is shown in Fig. [TOT a). This is quite promising: 
the N = 3 and 5 data coincide. 

Next, we can look at the top of the multiplet. Eq. [1] gives us a rescaled Lande interval 
rule: it says that the splitting between the J = N/2 and J = N/2 — 1 states is a constant, 
B, independent of J. Fig. [TOTb) shows this difference. It and panel (a) share the common 
N = 3 points, but the other points are different. The envelope of the curve is B(m q ). 

By design, these differences ignore the A term in Eq. [TJ To get it, we can look at the top 
of the multiplets, 

A = ^r M{N ' J = N/2) ~ ^r M{N ' J=N / 2 ~ x ) ( 22 ) 

or the bottom of the multiplets, 

A = ^M(N, J = N/2) - ±M{N, J = N/2-1). (23) 

Fig. [11] shows these two mass formulas. Again, they behave consistently. 

Figs. [10] and [TTJ show that the A and B coefficients in Eq. [TJ have typical hadronic sizes. 
Inserting a nominal lattice spacing, 1/a ~ 2100 MeV, we see that A = 400 MeV at small 
quark mass and is an increasing function of quark mass. In quark model language, A is 
the constituent quark mass, and its value and dependence on quark mass are both quite 
reasonable. B = 300 MeV and falls with energy. Again this is a typical hadronic scale. 

These values address an old conundrum of large-A phenomenology: the mass difference 
of the nucleon and A is supposed to be order Aq C d/N. It is measured to be about 300 
MeV, which is itself order Aqcd- The plots of mass differences show that the nucleon-A 
mass splitting is indeed well-described by large- iV QCD. Having data at several iV's as well 
as at several quark masses makes this result appear immediately. 
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FIG. 10: Exposing the B term of Eq. [1} (a) The J = N/2 vs J = N/2 — 1 mass difference in 
the SU(3), SU(5), and 577(7) multiplets, shown respectively as squares, diamonds, and octagons, 
(b) iV/3 times the J = 3/2 — 1/2 mass differences in the 577(3) and SU(5) multiplets, shown 
respectively as squares and diamonds. 
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FIG. 11: Exposing the A term of Eq. [TJ (a) The J = N/2 vs J = N/2 — 1 mass difference 
of Eq. [22] in the 577(3), 577(5), and 577(7) multiplets, shown respectively as squares, diamonds, 
and octagons, (b) Eq. [53] for the 577(3) and 577(5) multiplets, shown respectively as squares and 
diamonds. 
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SU(3) 517(5) 5*7(7) 



13 6.0175 

A = g 2 N 2.99 

r /a 5.36(6) 

n/a 3.90(3) 



17.5 34.9 
2.86 2.81 

5.20(5) 5.41(3) 

3.77(3) 3.91(2) 



rQ^fd 1.175(3) 1.172(3) 1.167(2) 
0.856(5) 0.850(4) 0.845(2) 



TABLE I: Bare parameters and observables from potentials. 



k am q amps amy amsiJ = 3/2) arag(J = 1/2) 



0.1230 0.119 0.525(2) 


0.617(3) 


0.996(6) 


0.942(5) 


0.1240 0.089 0.449(2) 


0.561(3) 


0.915(7) 


0.846(5) 


0.1250 0.055 0.362(2) 


0.505(4) 


0.835(8) 


0.744(6) 


0.1253 0.046 0.333(2) 


0.491(5) 


0.804(9) 


0.700(6) 


0.1260 0.029 0.255(3) 


0.457(6) 


0.747(12) 


0.613(7) 


0.1265 0.014 0.184(4) 0.427(10) 


0.703(13) 


0.516(12) 



TABLE II: Spectra from SU(3) simulations. 



V. CONCLUSIONS 



The qualitative features of large - A" QCD phenomenology are easy to observe. The shape 
of the potential (characterized by the scaling combination riy/a) is N - independent. Meson 
masses show little A" dependence - itlhTi is ^-independent. Vacuum - to -hadron matrix 
elements scale as y/~N. Two-flavor baryons show a prominent rotor spectrum. This seems 
to be true for both the bottom of the spectrum (low J) and the top (J = A" states). 

Naively, one might expect that the wave function of the baryon would be N - independent. 
Lattice simulations do not generally directly reveal wave function information, but one might 
expect that the same interpolating fields might behave the same for different iV's. That does 
not seem to be the case; larger A" seems to prefer larger Gaussian trial wave functions. But 
perhaps even A^ = 7 is not such large A" and one should push farther. 

Readers who do lattice simulations can undoubtedly list many obvious extensions of this 
project: smaller lattice spacing, to do an honest extrapolation to the continuum limit, bigger 
volumes to check that the answers are trustworthy, smaller quark masses, and bigger iV's 
are obvious technical improvements. Giving the two fermions different masses, or better 
yet, considering flavor S77(3), would allow more tests of large- A" spectroscopy. One might 
really want to test whether dynamical fermions become less important at large N. Writing 
the code to simulate SU(N) fundamental fermions is straightforward; running it might be 
costly, and seeing the expected small differences shrink as A" rises would be even more costly. 

The continuum literature of large - N baryons is more than thirty years old, and I clearly 
have only scratched its surface in this paper. Several simple ingredients were useful: having 
data at several iV's, at several J's for each N, and having data at many quark masses. Most 
of the continuum literature I have read restricts itself to statements about J = 1/2 and 3/2 
- presumably because that is all that exists in experimental data. Predictions for any J can 
challenge the lattice and would be candidates for future studies. 
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K arriq amps amy ame(J = 5/2) ame(J = 3/2) o)?ib(J = 1/2) 



0.1240 0.127 0.565(1) 


0.655(1) 


1.864(6) 


1.814(5) 


1.786(5) 


0.1250 0.099 0.488(1) 


0.593(1) 


1.708(6) 


1.650(6) 


1.616(5) 


0.1260 0.070 0.403(1) 


0.532(2) 


1.566(7) 


1.488(6) 


1.446(6) 


0.1265 0.055 0.356(1) 


0.500(2) 


1.486(7) 


1.404(6) 


1.356(6) 


0.1270 0.041 0.302(2) 


0.469(3) 


1.419(8) 


1.325(7) 


1.270(7) 


0.1275 0.026 0.240(2) 0.440(4) 


1.357(10) 


1.243(9) 


1.177(8) 



TABLE III: Spectra from SU (5) simulations. 



k am q amps amy amp(J = 7/2) a7ns(J = 5/2) ame(J = 3/2) 



0.1260 0.115 0.565(1) 


0.663(1) 


2.681(11) 


2.649(8) 


2.612(8) 


0.1270 0.088 0.488(1) 


0.603(1) 


2.483(9) 


2.420(10) 


2.386(9) 


0.1280 0.062 0.401(1) 


0.537(2) 


2.279(11) 


2.215(11) 


2.174(10) 


0.1290 0.036 0.299(1) 


0.471(3) 


2.082(15) 


1.994(12) 


1.929(11) 



TABLE IV: Spectra from SU(7) simulations. 
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